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CLOSED-LOOP SEPARATION CONTROL USING OSCILLATORY FLOW EXCITATION 


BRIAN G. ALLAN* * * § , JER-NAN JUANG*, DAVID L. RANEY*, AVI SEIFERT 8 , LATUNIA G. PACK', AND DONALD E. BROWNlI 

Abstract. Design and implementation of a digital feedback controller for a flow control experiment was performed. 
The experiment was conducted in a cryogenic pressurized wind tunnel on a generic separated configuration at a chord 
Reynolds number of 16 million and a Mach number of 0.25. The model simulates the upper surface of a 20% thick 
airfoil at zero angle-of-attack. A moderate favorable pressure gradient, up to 55% of the chord, is followed by a 
severe adverse pressure gradient which is relaxed towards the trailing edge. The turbulent separation bubble, behind 
the adverse pressure gradient, is then reduced by introducing oscillatory flow excitation just upstream of the point of 
flow separation. The degree of reduction hi the separation region can be controlled by the amplitude of the oscillatory 
excitation. A feedback controller was designed to track a given trajectory for the desired degree of flow reattachment 
and to improve the transient behavior of the flow system. Closed-loop experiments demonstrated that the feedback 
controller was able to track step input commands and improve the transient behavior of the open-loop response. 

Key words, active flow control, experimental, fluid mechanics, separation control, feedback, closed-loop 

AMS subject classification. Applied Numerical Mathematics, Controls 

1. Introduction. Experiments conducted by Seifert et al. [6, 7, 8] have shown that the introduction of oscillatory 
flow excitation into a separated flow region can fully reattach the flow or reduce the degree of flow separation. These 
experiments were performed on various airfoils at low and high Reynolds numbers. They demonstrated that the lift can 
be significantly enhanced when introducing periodic excitation upstream of the separated flow region. This enhance- 
ment is achieved by exciting the flow at frequencies that generate 2—4 spanwise coherent vortices over the length of 
the separated region. This periodic excitation increases the momentum transfer across the shear-layer, enhancing its 
resistance to flow separation due to curvature and adverse pressure gradients. The experiments conducted by Seifert 
and Pack [8, 9] also demonstrate that there is a range in which the response of the flow is proportional and even linear 
to the magnitude of the periodic excitation. These flow experiments have been conducted in an open-loop, or man in 
the loop, fashion. The goal of the present research is to add feedback control to the experimental setup used by Seifert 
and Pack [8], This controller will then vary the magnitude of the oscillatory excitation in order to achieve a desired 
degree of flow attachment characterized by a surface pressure gradient. 

2. Experimental Setup. The experiment was conducted in the 0.3-meter Transonic Cryogenic Wind Tunnel at 
NASA Langley Research Center. This tunnel is a fan driven, closed loop system, with a 0.33 m by 0.33 m test cross 
section and uses gaseous nitrogen (GNj) as the test medium. The tunnel operates at stagnation pressures ranging from 
1.2 bar up to 6 bar and total temperatures from 78 K up to 327 K [4, 5]. A fully automatic control system maintains 
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Fig. 2.1. Hump configuration. 


the test conditions, providing a high level of repeatability. 

2.1. Hump Model. The test model simulates the upper surface of a 20% thick airfoil that is a variation on the 
Glauert Glas II airfoil [8], The model contour generates a mild favorable pressure gradient until x/c = 0.55. This is 
followed by a severe adverse pressure gradient that is relaxed towards the trailing edge as shown in Fig. 2.1. Without 
control, the flow separates at the slope discontinuity, x/c « 0.66. Two alternative blowing slot locations are available: 
x/c = 0.59 and 0.64 (both shown on Fig. 2.1). The position of the upstream slot was selected so it would be upstream of 
the expected position of a shock wave. All of the experiments conducted in this investigation used the downstream slot 
at x/c = 0.64. The slots were about 0.25% chord wide (0.50 mm ± 20%), and allowed an ahuost tangential downstream 
introduction of momentum (the slots are inclined at 30° to the surface because of manufacturing considerations, see 
Fig. 2.1). 


2.2. Oscillatory Blowing Valve. A rotating, siren type, valve was used to generate pressure oscillations inside the 
model cavity. The oscillatory blowing valve was capable of generating frequencies up to 800 Hz, and for safety reasons 
was rated to 300 psi. GNo was supplied to the valve by converting a portion of the liquid nitrogen (LNi) available for 
operating the tunnel using an ambient temperature vaporizer. The use of ambient temperature GNj simplified the valve 
design. A pressure regulator was used to control the GNj entering the valve and the variable speed drive of the valve 
motor was used to control the frequency of the pressure oscillations. The oscillatory blowing valve was attached to 
the right tunnel plenum door at the center of rotation of the turntable (Fig. 2.2). The outlet of the oscillatory blowing 
valve was connected to the backside of the model cavity. The 76.2 mm inner diameter pipe, exiting the oscillatory 
blowing valve, was converted to a rectangular cross section and then split into five square diverging ducts that fit into 
the 330mm by 25.4mm entrance to the model cavity. Seven suction ports were positioned between and around the five 
square ducts at the exit of the manifold. These ports were connected to the left side tunnel boundary-layer removal 
system. The digital valves of the left side boundary-layer removal system were used to control the steady flow rate out 
of the manifold shown in Fig 2.2. This flow is driven by the tunnel static pressure which is higher than the ambient 
pressure. Any relevant combination of steady and oscillatory excitation could be generated with this type of control. 
Steady suction or steady blowing could be applied in a similar maimer by holding the oscillatory blowing valve in the 
fully open position and varying the inlet or the exhaust mass flow rates. 
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FIG. 2 . 2 . Experimental setup. 


2.3. Instrumentation. The model surface pressures were measured at 60 locations by the facility data acquisition 
system. Normal force, moment, and form-drag were calculated from these measurements. Turbine flow meters on the 
inlet of the oscillatory blowing valve and on the exhaust side of the model cavity (through the manifold) were used to 
measure the flow volume entering and exiting the model cavity. The volume flow rate measurements were combined 
with temperature and pressure measurements to determine the mass flux entering and exiting the model cavity (see Fig. 
2.2). The steady mass flux in or out of the slot could be determined by subtracting the two mass fluxes. A dynamic 
pressure transducer, flush-mounted in the model cavity, was used to measure the pressure fluctuations produced by 
the oscillatory blowing valve. A Thermocouple was installed inside the cavity to allow (using also the tunnel static 
pressure) calculation of the GNi density in the cavity. The normalized cavity pressure fluctuation, p'/p„, were shown 
to be directly related to the excitation unsteady momentum, <c„ >, that is the leading control authority of the periodic 
excitation when using optimal frequencies [9]. 

2.4. Closed-Loop System. The block diagram in Fig. 2.3 represents the setup of the experimental closed-loop 
system in the laboratory. This diagram shows the interconnection between the PC, actuator, and sensors. The amplitude 
of the control jet can be modulated using the PC which sends an IEEE pressure command, Pj, to the mensor unit. 
The mensor unit then outputs a pressure equal to Pj, to the dome pressure regulator, which supplies the GNj to the 
oscillatory blowing valve. Figure 2.3 shows how the digital exhaust valve is controlled by the PC which sends a second 
IEEE command, m ec , to the digital exhaust valve. This IEEE command is received by a unit which will open and close 
a desired combination of solenoid valves equalizing the mass flow rate entering and exiting the cavity of the model. 

The pressure gradient, dCp /dx, was measured using a differential pressure transducer. This sensor measured the 
pressure difference between stations, x/c = 0.952 and 0.489. This pressure gradient is used to characterize the degree 
of flow reattachment on the hump model and will be referred to as the pressure recovery parameter. The output voltage 
from the transducers are sent to a low pass filter, amplified, and then sampled by the DSP board on the PC. The signal 
from the dynamic pressure sensor, inside the cavity of the hump model, is passed through a RMS to DC converter 
and then amplified before being sampled by the DSP board. The DSP board samples these two signals at a rate of 
100 Hz. This data can either be uploaded as raw sampled data to the PC or averaged over a desired length of tune and 
then uploaded to the PC. In the closed-loop experiments conducted here, the sampled data was averaged over a 0.5 s 
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Fig. 2.3. Block diagram of experimental system. 


interval and then uploaded at every second from tire DSP board to the PC. 

3. System Model. A linear system model was developed for the feedback control design process. The actual 
system has some nonlinear effects which are neglected in the linear system model but will be incorporated in a non- 
linear model. The nonlinear model will be used to evaluate the feedback controller developed from the linear system 
model. 


3.1. Continuous System Model. The dynamics of the flow system were modeled by lumping the system com- 
ponents into two separate systems as shown in Fig. 3.1. The first system models the dynamic response from the IEEE 
command, Pj—P s , to the RMS pressure fluctuations in the cavity, p' / p TO . This system will represent the dynamics 
of the actuator system. The second system represents the dynamic response of the flow system and models tire re- 
sponse of the pressure recovery parameter, dCp /dx, to the cavity RMS pressure fluctuations, ///p„. Both systems are 
modeled using a simple second-order damped system. The second-order model was chosen since the dynamics of the 
experimental system are unknown and the dynamics of most mechanical systems can be approximated by this type 
of general model. If the second-order model is determined to be insufficient a model of increased complexity can be 
used. 

The dynamics of the actuator system are modeled by the simple second-order system 


(3.1) 


£_ 

dt 2 


P'(t) 


— 2Ci w «i . 
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where is the damping ratio and <n„i the natural undamped frequency. The constants m.\ and b\ are the slope 
and intercept of the linear fit to the steady-state data of p' / p TC as a function of Pj—P s . The steady-state relation 
{p'/p°o)ss = f{[Pi—Ps]ss), is then approximated by the linear fit 

(3.2) ) =m i (P i -P s )„ + bi 

\ P°° / ss 


where (p 1 /p°o) S s and (P,—P s )ss are the steady-state values of (p 1 / p TC ) and (P, — P s ) respectively. Similarly, the dynamic 
response of the pressure recovery parameter, dCp / dx, to the cavity pressure fluctuations, p'(t ) /p M , is modeled by the 
second-order system 


(3.3) 
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System Model 


u = Pi 



Fig. 3.1. Two system model 


The constants mo and b 2 are the slope and intercept of the linear fit to the steady-state data for dCp /dx as a function 
of p' / Poo. The steady-state function (Cp) ss = g([p' / p^\ss) is then approximated by the linear fit 


(3.4) 


0 dC P /dx) ss = m 2 (p'/poo)ss+b 2 


Next, define the state vector x = [x 4 ,x 2 ,X 3 ,X 4 ] T = [p 1 /p°°, d /dl(p' / p^), dCp/dx , d /dt(dCp/ dx)] T . Using this state 
vector, Eq. 3.1 and Eq. 3.3 can be expressed as 
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Substituting p'/poo = xi for the input term results hi 


(3.6) 
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Defining u = Pj—P s , the equation above can be expressed in the matrix form 

(3.7) x = Ax + Bu + f 

where the vector f is a drift (or bias) term. The output vector y = [y i ,y 2 ] J is defined as 


(3.8) 


y = 


p'/poo 


o 

o 

o 

1 

dCp/dx 


0 0 1 o J 


x = Cx 


The drift term f can be removed by considering the steady-state solution to Eq. 3.7 

(3.9) 0 = Axjj + f 

where u ss and x ss are the steady-state input and state vectors respectively. A new state variable x is now defined as 

(3.10) x = x-x ss 


Substituting this relation into Eq. 3.7 and using Eq. 3.9 results in 
(3.11) ic = Ax + Bv 
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FIG. 3.2. Simplified block model 


where v = u — u ss . Using the definition of x, the output equation for y becomes 


(3.12) 


y = Cx + Cx ss 
= y + Cx ss 


where y = Cx. The new system without the drift tenn can be summarized as 


(3.13) 

where 


x = Ax + Bv 


y = Cx 


x = x+x ss 

(3.14) y = y + Cx ss 

U=V+ Uss 

The system in Eq. 3.13 is expressed as the simplified block, P, shown in Fig. 3.2 where P has one input v and two 
outputs yi and V 2 - This linear system without the drift tenn can now be used in the feedback control design. 

3.2. Discrete System Model. Sampling of the analog voltages from the pressure transducers was done digitally, 
making our control system discrete. In order to account for the discrete sampling in our controller, a discrete control 
design was performed. To design our controller in the discrete domain, the continuous system given in Eq. 3.13 was 
converted to a discrete system using a zero-order hold. 


Xd (k + 1 ) = A d x d ( k ) + B d \ d ( k ) 

yd(£) = Qx d (£) 


The discrete state variable is defined as, Xd(k) = x(kT), where T is the sampling period. Using a zero-order hold on 
the input, the matrices A d and B d are given by 


(3.16) 


A d = e AT 

B d = (A d -I)A~ l B 


where C d —C and A is assumed invertible [3]. 

The parameters for the A and B matrices were determined using steady-state and dynamic open-loop experimental 
data. The experimental data used to determine the steady-state parameters are shown in Figs. 3.3 and 3.4. Figure 3.3 
shows the experimental data for the steady-state relationship between the IEEE command P,. minus the static wind 
tunnel pressure P s , to the cavity pressure fluctuations p' divided by the free stream density p TC . Tins figure shows that 
the linear fit represents the experimental data reasonably well, for P,—P s above 5. Below P,—P s = 5. the cavity pressure 
fluctuations do not respond to the low P, input. Notice that some of the experimental data points do not fall on the 
lmear fit. This will result in steady-state errors in our linear model near these points. Fig. 3.4 shows the experimental 
data for the steady-state relation between the cavity pressure fluctuations and the pressure recovery parameter, dCp / <lx. 
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FIG. 3.3. Steady-state effect of Pj—P s on p* / p c 



FIG. 3.4. Steady-state effect, of p' /p^ on dCp/dx. 


Table 3.1 

Dynamic parameters for system model. 


tn\ 

55.77 

bi 

-288.2 

(Dtfi 

0.5000 

111 

0.3400 


m 2 

4.956- 10“ 4 

b 2 

1.781 

t0„ 2 

80.00 

*|2 

1.000 


This figure shows a change in the gam for values of p' / < 400. The linear fit is made to the data above p' / p TO = 400 
as tire controller will typically be operating at cavity pressures fluctuations in this range. The steady-state gams for file 
linear system are given by the parameters m 2 , and b 2 shown in Table 3.1. 

Figure 3.5 shows the hump model pressure distribution with and without control. The baseline pressure distribu- 
tion, given by the P, —P s = 0 line, shows a large separation bubble downstream of the slot. For file baseline pressure 
distribution, the pressure recovery parameter, dCp/dx, is equal to 0.77. Applying periodic excitation (Pj — P s > 0) 
decreases the size of the separation bubble, as shown by the pressure distribution in Fig. 3.5. The P, — P s = 10 line 
shows an increase in the pressure recovery parameter, dCp/dx, of 1.6, as compared to 0.77 given by the baseline 
pressure distribution. Increasing P, — P s to 20, produces a dCp / dx value of 1 .85, indicating a further size reduction of 
the separation bubble. 

Figures 3.6 through 3.9 show the open-loop responses of the experimental system to four different P/—P s step 
input commands. These open- loop responses were used to determine the dynamic parameters in Table 3.1 so that 
the linear model response would match the experimental results. By changing the values of to,, 1 , £1, ra„2, and 'Cp_ 
the dynamic characteristics of the linear system model could be fitted to the data. The open-loop experimental data 
showed that the actuator system was much slower than the response of the flow system. The frequency response of the 
flow was fixed at oo„2 = 80 with a damping ratio of Lp_ = 1. The parameters to„ 1 and were then adjusted to fit file 
experimental data. The oscillations in the experimental data for dCp/dx were not modeled by the linear model. These 
small oscillations are believed to be produced by small fluctuations in the wind tunnel test conditions. 

Figure 3.6 shows a comparison of the experimental data to file linear model response to a step input from Pj—P s = 
10 to 15. Here the linear model over predicts the initial steady-state cavity pressure fluctuations and pressure recovery 
parameter for P,—P s = 10 and under predicts the cavity pressure fluctuations at the final steady-state at Pj—P s = 15. 
This steady-state offset is a result of the linear fit to the steady-state data. Figure 3.6 does show a good fit to the 
dynamic response of the experimental system for the cavity pressure fluctuations. The comparison for the pressure 
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FIG. 3.5. Steady-state surface pressure profiles for Pj — P s = 0, 5, 10, and 20 


recovery parameter shows good agreement for the transient and the mean steady-state values. Figure 3.7 shows a good 
prediction of the open-loop dynamics but under predicts the cavity pressure fluctuations. This under prediction by the 
model is again due to the lineal' fit of the steady-state experimental data. Figure 3.3 shows that the two experimental 
points at Pj—P s = 15 and 17 are slightly higher than the linear fit. This difference then shows up in the open-loop 
simulations as an offset. Figure 3.8 shows a slight over prediction of the steady-state pressure recovery parameter but 
does capture the dynamics of the step input. Figure 3.9 shows that the model has an offset in both the cavity pressure 
fluctuations and pressure recovery parameter. This offset is due to an error in the steady-state prediction as can be seen 
in the linear fit to the steady-state data in Figs. 3.3 and 3.4 for Pj—P s = 25. Overall the linear model does a fair job 
of predicting the dynamic response of the experimental system and has a slight offset in the steady-state gams of the 
system. 

Using file parameters in Table 3. 1 the matrices A and B, for the continuous system in Eq. 3.6, are given as 


(3.17) 



0 

1.0 

0 

0 


0 


-0.25 

-0.34 

0 

0 


13.943 

A = 

0 

0 

0 

1 

B = 
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1.798 

0 

-6400 

-160 


0 


Using a sample period of T = 1 .O.v, the matrices A,/ and B,j for the discrete system, using Eq. 3.15, are 


(3.18) 


Ad = 


0.89030 

0.81292 

0 

0 ’ 


6.1180 

-0.20323 

0.61391 

0 

0 


11.334 

2.5158 - 10“ 4 

2.2401- 10" 4 

0 

0 

Bj = 

1.6407 -10- 3 

—5.6003 • 10 -5 

1.7536 - 10" 4 

0 

0 


3.1240- 10~ 3 


Note that since the sampling frequency is much lower than the frequency of the flow system, the A,/ matrix results 
in the last two columns being zero. Since the sampling frequency is so low these states could be excluded but been 
included for completeness. Future actuators will have an unproved response in the high frequency range where file 
inclusion of the flow dynamics will be necessary. 

4. Feedback Control Design. A feedback controller was designed for this system using a simple proportional- 
integral-derivative (PID) controller. The advantage of using a PID controller is that it can be tuned to give good 
performance results where only the dominate system tune constants are known [2], 

Figure 4. 1 shows a block diagram of a basic PID feedback controller for a continuous closed-loop system. This 
figure shows the output, > (/), being subtracted from the desired output. y f /„((), to produce an error signal. e(t). Given 
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the error signal, e(t), the input, u(t), can be computed by 

(4. 1) w(0 = K P e(t) + K, £ e(t)dt+K D ^ 

The constant Kp can be adjusted to vary the amount of proportional feedback for the controller. Likewise, the constants 
Kj and Kp, are used to vary the degree of integral gain and derivative gam respectively. 

Increasing the proportional gain will typically reduce the steady-state error and settling time of the closed-loop 
system. However, large Kp values usually produces instability hi the closed-loop system. For most systems, there is 
an upper limit on the gam Kp in order to achieve a well-damped, stable response. This upper limit on Kp may still 
result in unacceptable steady-state errors, putting a limit on how much the steady-state errors can be reduced using 
proportional feedback alone. 

To improve the steady-state accuracy, without needing a large proportional gain, integral control can be used. 
Integral control will reduce or eliminate steady-state errors at the cost of reduced transient response performance. 
When the steady-state error is small, the proportional gain times the error will produce a small control input, assuming 
that the proportional gam can not be arbitrarily large. This small input from proportional feedback alone, will typically 
be too small to change the steady-state error. Integral gain on the other hand, will integrate the small error over time, 
producing an increasing control input until the system responds. This input from the integral gam will eventually 
reduce or eliminate the steady-state error. Like the proportional gam, the integral gain, Ki, also has an upper limit in 
order to achieve a well-damped, stable response. 

Derivative control can be added to increase the damping and increase the stability of the closed-loop system. Since 
the derivative of the error represents the slope of the error, <?(/ ), the derivative feedback is essentially an anticipatory 
type of control. Normally, if the slope of the error is large, an overshoot will subsequently occur. The derivative 
control can be used to predict an overshoot by measuring the slope of the error and making a proper correcting effort 
before the overshoot actually occurs. 


4.1. Discrete PID Controller. The PID controller was implemented using a computer which resulted in the need 
for a digital controller. To account for the discrete nature of the digital controller in the control design, the controller 
was designed in the discrete domain. 

Approximation of the continuous integral controller was done using a trapezoidal method. Let gi ( kT ) be the 
numerical approximation to the integral of the error at t = kT, where g/(0) = 0. Using the trapezoidal method, gi(kT) 
can be written as 

(4.2) gI (kT) = gI [(* - 1 )T] + | [e(kT) + e((k - 1 )T)] 


Taking the z-transfonn on both sides of Eq. 4.2 results in the transfer function 


(4.3) 


G/(-) 


E(z 


n=+ 1 ) 

2 ( 2 - 1 ) 


where M/ (z) is a discrete transfer function from the error to the integral of error. The discrete approximation to the 
derivative was made using a backward difference method. Let go(k'i) be the approximation to the derivative of the 
error at t = kT resulting in 

(4.4) gD (kT) = ^{e(kT)-e[(k-l)T\} 


Taking the z-transform on both sides of Eq. 4.4 results in the transfer function 


(4.5) 


M d (z) 


Gd{=) 
E( z) 
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where Mo ( 2 ) is a discrete transfer function from the error, to the derivative of the error. Combining the proportional, 
integral, and derivative controllers, results in the transfer function from the error, E(z), to the control input U (z) 
becomes 


(4.6) H(z) = tttt = Kp +K,M,(z) + K d M d (z) 

t(z) 

where H ( 2 ) is the discrete transfer function for the PID controller from the error signal to the control input. Substituting 
Eq. 4.3 and 4.5 and simplifying results in 


(4.7) 


H(z) 


K, + ^ + ±W + 


—Kp + 


IK, 


z{z- 1) 


2Ko\ Ko 
T / + T 


Figure 4.2 shows a block diagram of the transfer function in Eq. 4.7. 


4.2. Discrete PID Control Design. The block diagram in Fig. 4.3 shows the closed-loop system with the discrete 
PID controller, 7 /( 2 ), and the linear system model, P(z). The closed-loop system also includes an external disturbance, 
w, which can be used to model the effects of the pressure oscillations created by small fluctuations in the wind tunnel 
test conditions. The effects of these fluctuations can be seen in the open-loop experimental data in Figs. 3.6 through 
3.9. The objectives of the control design is to track a desired dC P /dx command while minimizing the overshoot. Since 
the actuators for this system were fairly slow, minimizing the effects of wind tunnel fluctuations was not possible and 
therefore not a control objective. The response time of the controller was considered a secondary objective and was 
relaxed in favor of smoother transient behavior. 

Figure 4.4 shows a root locus plot in the z-plane for the poles of the discrete closed-loop system. The plot also 
shows the lines of constant damping, £, and natural frequency, o)„ , transformed from the s-plane to the z-plane. Poles 
inside of the unit circle, in the z-pane, are stable and poles outside are unstable. The open-loop system has four poles 
which are shown in Fig. 4.4 as open circles. The two high frequency poles in the s-plane are mapped to the origin, 
since the sample frequency is well below the frequency of these poles. Increasing K P while holding K, = Ko = 0, 
moves the two low frequency open-loop poles to a higher natural frequency, while decreasing the damping ratio. The 
closed-loop poles cross the unit circle at Kp = 117 where lar ger values of Kp result in an unstable closed-loop system. 
Increasing the integral gain, K,, while holding Kp = Kp = 0, also moves the low frequency poles towards the unit 
circle. At K, = 12.7, the poles cross the unit circle where larger values of K ) make the closed-loop system unstable. 
The integral controller also introduces a fifth pole which moves along the real axis, from one to zero, as K, increases. 

The time response of the closed-loop system, for increasing values of Kp and Kj, ar e shown in Figs. 4.5 through 
4.8. The figures show the response of the closed-loop system, where the controller is turned on at t. = 0 with a desired 
dCp /dx = 1.8. Figures 4.5 and 4.6 show how using only proportional gain results in a significant steady-state error. 
This error is a result of the constant control input needed in order for the system to reach the desired steady-state point. 
Since the proportional gam is a product of the error, the error could never go to zero while generating the constant 
control signal needed for the new steady-state point. By increasing the gain Kp, the error is reduced only slightly at the 
cost of reduced damping. This reduced damping effect can be explained by the root locus plot in Fig. 4.4 where the 
closed-loop pole is moved closer to the unit circle as Kp is increased. Figures 4.7 and 4.8 show the effect of only using 
an integral controller for a step input command. These figures show how the integral controller dramatically reduces 
the steady-state error as compared to the Ko gam. Unlike the proportional gain, the integral gam is able to generate a 
constant input signal in order to reach the new steady-state point. By increasing K,, the response time decreases but 
at the cost of decreased damping. Like the proportional controller, increasing K, moves the two poles toward the unit 
circle decreasing the damping for these poles. 
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To improve the damping of the closed-loop system for high gains, derivative feedback gain can be added. However 
this gam can cause problems when there is significant noise in the system as is the case here. The noise in the system 
was shown in the open-loop experimental data in Figs. 3.6 through 3.9. As a result of the noise content in our system, 
the derivative feedback gam was not used. However if minimizing the response time was the main control objective 
then derivative feedback might be needed. 

To see how the Kp and Kp gams are effected by the wind tunnel fluctuations, a Bode plot for the transfer function 
between the disturbance input, w, and error, e, are plotted in Figs. 4.9 and 4.10. Figure 4.9 shows how increasing 
the K[ amplifies noise signals in the 0.03 to 0. 1 Hz range and remains fairly flat above 0. 1 Hz. This is good since the 
large amplitude oscillations due to the wind tunnel noise fall in the 0.17 to 0.12 Hz range. Figure 4.10 shows that 
increasing Kp will increase the error signal for noise in the 0.08 to 0.3 Hz range. This means that using a proportional 
gain will amplify the oscillations produced by the fluctuations in the wind tunnel test conditions. Therefore in order 
not to amplify these oscillations, proportional gam was not used. 

5. Nonlinear Simulation. A nonlinear model of the experimental system was developed using the SlMULINK 
toolbox which operates in the Matlab computing environment [1], This model includes the nonlinear effects of the 
steady-state gams as shown in Figs. 3.3 and 3.4, as well as the limiters applied to the control input. This simulation also 
accounted for the discrete time sampling and averaging of the pressure sensor signals by the DSP board. Figure 5.1 
shows file mam block diagram of the experimental system which simulates file closed-loop step response of dCp / dx to 
a step input command in file desired dCp /dx value. This simulation model was used to evaluate the feedback controller 
after it was designed using the linear system model. 

The discrete PID controller is shown as a single block in Fig. 5.1 with two inputs and one output. The contents of 
this subblock are shown in Fig. 5.2 where the two inputs are subtracted to make an error signal which is then passed to 
the proportional, integral, and derivative gains. These gains are summed up to make file control output signal Pj—P s . 
The Pj—P s signal is then passed through the discrete rate and saturation limiters shown in Fig. 5.3. These limiters were 
used to keep the flow system in a safe operating range. The output from the limiters is then sent to the flow system 
model shown in Fig. 5.4. This system models the nonlinear response of the actuator from the input signal Pj—P s , to 
the output signal dCp /dx. A lookup table was used for the nonlinear steady-state gains and the dynamics of file system 
were modeled using the second-order system described in section 3. 

The closed-loop model shown in Fig. 5.1 uses the sensor measurement of the pressure recovery parameter, 
dCp /dx, which is sampled by the DSP board. The DSP board samples voltage of the differential pressure trans- 
ducer at a frequency of 100 Hz. This signal is sampled and averaged over a tune period of 0.5 s. The averaged signal 
is then sent to the discrete controller at a rate of 1 Hz. The averaging procedure is modeled by file averaging block in 
Fig. 5.1, winch is shown in detail in Fig. 5.5. The system hi Fig. 5.5 starts summing up the 100 Hz sampled data when 
given a trigger command. The trigger in Fig. 5.1 is a square wave which is on for 0.5 5 and then off for 0.5 s. 

Figures 5.6 through 5.9 show a comparison of four different open-loop step responses between the experiment, 
linear system, and the nonlinear model. The nonlinear model shows a better prediction of the steady-state gams since 
it is able to interpolate from the experimental data. An overall comparison of the dynamic response of the systems 
show that the linear model and the nonlinear model are nearly identical as might be expected. 

6. Results. Based on the PID discrete control design, it was predicted that using only integral control would best 
meet the design objectives for the feedback controller. The linear simulations using the integral controller in Fig. 4.8, 
showed the closed-loop response for K/ = 5,6, and 7. Setting the integral gain K[ = 6.946, the feedback controller 
was then implemented on the experimental closed-loop system. The results from this experimental data are shown in 
Fig. 6.1. The test was performed by setting the value of P, —P s = 10 and then turning on the controller at t = 0 given 
a desired pressure recovery parameter, dCp/dxj es = 1.9. The experimental data shows that the controller performed 
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well and had very little overshoot and a good response time. Figure 6. 1 also shows the predicted response of the linear 
and nonlinear models. Both models were able to predict the experimental closed-loop response relatively well and the 
transient response for dCp/dx was predicted very well by the linear model. The prediction by the nonlinear model 
for the dCp / dx transient was good, but had a slightly larger overshoot than the experimental results. Both the linear 
and nonlinear models showed larger oscillation in the cavity pressure fluctuations p' / and inlet pressure P, —P s as 
compared to the experimental data. Figure 6.2 shows a comparison between the open-loop response and closed-loop 
response for the experimental system. The open-loop response was performed by using a step input in Pj — P s using 
the initial and final values given by the closed-loop simulation. This figure shows that dCp/dx and p' / for the 
open-loop response, has a faster response tune but at the cost of a larger overshoot. This figure also shows that there 
is a slow transient in file flow system. At t = 25 s, dCp/dx has reached its desired value of 1.9, yet the cavity pressure 
fluctuations are still being increased by file integral controller. On file other hand, the open-loop step input shows that 
dCp /dx is slightly higher than the desired value at I = 25 ,v and then slowly drifts to the desired value. This shows die 
advantage of the feedback controller since it is able to hold the desired pressure gradient while the flow system goes 
through a slow transient. 

Using file same integral gain, file controller was tested using a step down in file desired dCp /dx command. In 
this regime the steady-state gams start to deviate from the linear model fit as shown in Fig. 3.4. Figure 6.3 shows the 
closed-loop experimental data where the initial inlet pressure was set to, P, — P s =20 and the controller was turned 
on at t = 0 with dCp /dx = 1.64. The response of the closed-loop system displayed large oscillations in the dCp /dx 
signal, which settled down near t = 50 s. The nonlinear model also shows large oscillations which did not settle down 
where the linear model had mild oscillations for the closed-loop simulation. Neither model was able to predict file 
closed-loop response in this regime but they did predict the unsteady behavior for Kj = 6.946. The linear model is 
certainly not valid in this regime since it does not account for the steady-state gain shown in Figs. 3.3 and 3.4. The 
nonlinear model does include the nonlinear steady-state gain but is still unable to accurately model the closed-loop 
system. To improve the nonlinear model in this regime, more steady-state and dynamic data is needed. 

By reducing the integral gain, the closed-loop response in the nonlinear regime can be improved. Figure 6.4 
shows the closed-loop experimental data for K[ = 4.63 along with the linear and nonlinear simulations. The dCp /dx 
transient has improved but there ate still some mild oscillations in the experimental closed-loop system. The nonlinear 
model still predicts a continuing oscillation for this gain and the linear model is able to capture the average transient 
response but does not capture the oscillations seen in the experimental data. Reducing the integral gam to Kj = 2.315, 
the transient response is much smoother at the cost of a slower response tune as shown in Fig. 6.5. Both the linear 
and nonlinear models show good agreement with the experimental data. Figure 6.6 shows a comparison between 
the closed-loop and open-loop step input responses. This comparison shows how the open-loop response has a large 
overshoot in the cavity pressure fluctuations and some oscillations in the value of dCp / dx. The integral controller 
shows how it increases the damping effect of the transient response by shaping the input command P, — P s . This 
smoother inlet pressure command does not excite the higher frequency dynamics of the system like the open-loop 
step input. Figure 6.7 shows the effect of decreasing the integral gain for the experimental closed-loop system. This 
comparison shows how a smaller integral gain of Kj = 2.315 minimizes the overshoot and oscillations produced by 
the controller. 

7. Conclusions. This investigation demonstrated the incorporation of digital feedback control into an open-loop 
flow control experiment. The objective of the feedback controller is to hack a desired pressure gradient measured 
in the controlled flow region. This pressure gradient characterizes the degree of flow separation (or attachment) and 
hence the lift, drag, and moment forces acting on the generic hump model. The oscillatory excitation, which controlled 
the flow in the separated region, operated at 385 Hz. The degree of flow separation was then controlled by changing 
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the magnitude of the oscillatory excitation. The actuator system which changed the magnitude of the oscillatory 
excitation, had a frequency response of 0.05 Hz, which is much slower than the flow response. Therefore from the 
flow physics point of view, the changes hi the magnitude of the oscillatory excitation are performed in a quasi-steady 
manor. The dynamics of the open-loop system are then the dynamics of the actuator system and do not include 
dynamics of the controlled flow region. Experimental data from open-loop step inputs showed that the actuator system 
could be modeled by a simple second-order model. 

A PID control design approach was used since it could be tuned to give good performance results where only 
the dominant system tune constants are known. The PID control design showed that using only integral gam was 
the most effective as a result of the large time constant of the actuator system. The control design showed that 
adding a proportional or derivative gain would only amplify the oscillations produced by the small fluctuations in 
the wind tunnel test conditions. The closed-loop experiments demonstrated that the integral feedback controller was 
able to improve the transient response of the open-loop system by minimizing the overshoot seen in the open-loop 
experiments. Overall, the integral controller did a good job in tracking the desired pressure gradient in the controlled 
flow region. It should be noted that while the PID control design was reasonably simple, the application of a digital 
feedback controller to a flow control experiment is unique. This closed-loop experiment is the first in a series of 
experiments which will incorporate digital feedback control to existing NASA Langley flow control experiments. The 
knowledge and insights gained from this experiment will be valuable for future closed-loop flow control experiments. 
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FIG. 3.6. A comparison of an open- loop step input ofPj—P s 
from 10 to 15 for the experimental system and the linear discrete 
model 




FIG. 3.7. A comparison of an open-loop step input ofPj—P s 
from 15 to 11 for the experimental system and the linear discrete 
model 




Fig. 3.8. A comparison of an open-loop step input of Pi— P s 
from 17 to 20 for the experimental system and the linear discrete 
model 




Fig. 3.9. A comparison of an open-loop step input ofPj—P s 
from 20 to 25 for the experimental system and the linear discrete 
model 
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Fig. 4.1. Block diagram for PID controller. 



Fig. 4.2. Block diagram for digital PID controller. 
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Fig. 4.3. Block diagram of the closed-loop system. 
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Fig. 4.4. Root locus of closed-loop system in the z-ptane using proportional and integral feedback gain. The Kp lines show how the open-loop 
poles ( circles ) move using only proportional feedback. Similarly the Kj tines show how the open-loop poles move when using only integral feedback. 
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FIG. 4.5. Closed-loop simulation using the linear state model for FIG. 4.6. Closed-loop simulation using the linear state mode! for 

feedback gains Kp = 10, 20, 30 where K] = Kp, = 0. feedback gains Kp = 40, 50, 60 where K\ = Kp = 0. 
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Fig. 4.7. Closed-loop simulation using the linear state model for 
feedback gains Kj =2,3,4 where Kp =Kp = 0. 


FIG. 4.8. Closed-loop simulation using the linear state model for 
feedback gains Kj = 5,6,7 where Kp = Kp, = 0. 
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FIG. 4.9. Bode plot of closed-loop transfer function 
from the noise input w to the error e,for increasing \alues of 


FIG. 4.10. Bode plot of closed-loop transfer function 
from the noise input w to the error e,for increasing values of 

















FlG. 5.2. Block diagram of the PID controller. 
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FlG. 5.3. Block diagram of the rate and saturation limiters. 
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FIG. 5.4. Block diagram of the flow system model. 
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Fig. 5.5. Block diagram of the averaging block. 
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Fig. 5.6. A comparison of an open-loop step input ofP,—P s 
from 10 to 15 for the experimental system, the linear discrete model, 
and the nonlinear model 




FIG. 5.8. A comparison of an open-loop step input ofPj—P s 
from 17 to 20 for the experimental system, the linear discrete model, 
and the nonlinear model 




Fig. 5.7. A comparison of an open-loop step input ofP,—P s 
from 15 to 17 for the experimental system, the linear discrete model, 
and the nonlinear model 




Fig. 5.9. A comparison of an open-loop step input of P, — P, 
from 20 to 25 for the experimental system, the linear discrete model, 
and the nonlinear model 
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Fig. 6.1. Closed-loop response of the nonlinear model and 
the experimental system for (dCp / dx)j es = 1.9 where the controller 
was turned on at t = 0. The PID controller had the gains, Kp = 
0, K, =6.946, andK D = 0. 

Fig. 6,2. A comparison of the open-loop step input response 
to the closed-loop system using Kj = 6.946. 

OR 


PRr i 


on, 

Nonlinear Model 

r) 

Nonlinear Model 

^0' 

to 

n 

\ Experiment 

^ y ..................... 

— Linear Model 

— Experiment 


Time, sec 


Time, sec 



FIG. 6.3. Closed-loop response of the nonlinear model and the 
experimental system for (dCp fdx) des = 1.64 where the controller 
was turned on at t = 0. The FID controller had the gains, Kp = 
O.AT/ =6.946, andKo = 0. 


FIG. 6.4. Closed-loop response of the nonlinear model and 
the experimental system for ( dCp/dx)j es = 1.9 where the controller 
was turned on at t = 0. The FID controller had the gains, Kp = 
0,Kj = 4.64, andK D = 0. 














Fig. 6.5. Closed-loop response of the nonlinear model and 
the experimental system for (dCp / dx)j es = 1.9 where the controller 
was turned on at t = 0. The PID controller had the gains, Kp = 
0, K, =2.315, and K d =0. 
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Fig. 6.6. A comparison of the open-loop step input response 
to the closed-loop system using K / = 2.315. 



Fig. 6.7. A comparison of the experimental closed-loop transient responses for three K] values. 
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